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In the framework of a Varying Speed of Light theory, we study the eigenvalues associated with the 
Wheeler-DeWitt equation representing the vacuum expectation values associated with the cosmolog¬ 
ical constant. We find that the Wheeler-DeWitt equation for the Friedmann-Lemaitre-Robertson- 
Walker metric is completely equivalent to a Sturm-Liouville problem provided that the related 
eigenvalue and the cosmological constant be identified. The explicit calculation is performed with 
the help of a variational procedure with trial wave functionals related to the Bessel function of the 
second kind K v ( x ). We find the existence of a family of eigenvalues associated to a negative power 
of the scale. Furthermore, we show that at the inflationary scale such a family of eigenvalues does 
not appear. 


I. INTRODUCTION 

What is the Cosmological Constant? How can be computed? These are some of the many puzzling questions 
which are still unsolved. Basically the Cosmological Constant can be connected to the energy of the vacuum. How¬ 
ever, the absence of a complete Quantum Gravitational theory increases the number of questions instead of giving 
answers. General Relativity (GR) is the best theory explaining the behavior of the gravitational field including also 
the cosmological constant. However GR fails to describe the gravitational field in the quantum range. Despite of 
this problem, in GR there exists a quantization procedure known as the Wheeler-De Witt equation (WDW) [I] which 
encodes some aspects of the quantum properties of the gravitational field included the cosmological constant. We say 
“some”, because a complete solution of the WDW equation does not exist and one needs to reduce the degree of the 
difficulty by fixing a background and freezing some degrees of freedom. The WDW equation is the quantum version 
of the classical Hamiltonian constraint representing the invariance under time reparametrization. Its derivation is a 
consequence of the Arnowitt-Deser-Misner (ADM) decomposition (2] of space time based on the following line element 

ds 2 = 9uv ( x ) dx^dx v = (—iV 2 + NiN 1 ) dt 2 + 2Njdtdx J + gijdx' L dx\ (1) 

where N is the lapse function and N, t the shift function. In terms of the ADM variables, the four dimensional scalar 

curvature R can be decomposed in the following way 

U = R + K i:j K ij - (K) 2 - 2V„ (Ku" + a M ), (2) 

where 

K%j = ~2Jg \Pt9ij ~ N Mj — Ajp] (3) 

is the second fundamental form, K = g lJ K L] is its trace, R is the three dimensional scalar curvature and yfg is the 

three dimensional determinant of the metric. The last term in (0 represents the boundary terms contribution where 
the four-velocity vT is the timelike unit vector normal to the spacelike hypersurfaces (t=constant) denoted by £ t and 
aT = u a X/aU 11 is the acceleration of the timelike normal Thus 

C [N, Ni, 9ij ] = (72. - 2A) = [ KijK U - K 2 + R - 2A - 2V„ (Ku» + a")] (4) 
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represents the gravitational Lagrangian density where k = 87 tG with G the Newton’s constant and for the sake of 
generality we have also included a cosmological constant A. After a Legendre transformation, the WDW equation 
simply becomes 


WS> = 


(2k) G~ kl 


# R - 2A > 


'F = 0, 


( 5 ) 


where Giju is the super-metric and where the conjugate super-momentum it 13 is defined as 


SC 


S ( d t gi 


= (g i3 K - K ij ) 
yy ’ 2 k 


( 6 ) 


Note that H = 0, represents the classical constraint which guarantees the invariance under time reparametrization. 
The other classical constraint represents the invariance by spatial diffeomorphism and it is described by n^ 3 = 0, 
where the vertical stroke “|” denotes the covariant derivative with respect to the 3 D metric (jij . Solving Eq.© allows 
to extract information on the early universe and on the cosmological constant. Of course, the form of the solution 
is depending on the background one considers. In this paper, we fix our attention on the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric without matter fields. To the reader, this choice could seem a restriction, however 
one has to think that in the very early universe, before the inflationary phase, it is likely that all the quantum 
information can be carried on by the gravitational field, because of its non-linear nature. However, even in this 
simplified vision many problems arise, especially for the inflationary epoch. In recent years, the idea of modifying 
GR to cure some of its diseases has been considered. From one side, the so-called f (R) theories have been taken 
under examination to cure some problems in the infrared (IR) region[3j] and on the other hand modifications on short 
scales allowing a power-counting ultraviolet (UV) renormalizable have been proposed by Horava motivated by the 
Lifshitz theory in solid state physics [UfE]. This theory is dubbed as Horava-Lifshitz (HL) theory and should recover 
general relativity in the IR limit. Nevertheless the price to pay to obtain a renormalizable theory in the HL proposal 
is that we have no general covariance or, in other words Lorentz symmetry is broken. Another proposal which distorts 
gravity in the UV is Gravity’s Rainbowjjj]. Gravity’s Rainbow has some appealing features to explain inflation^. 
In a series of papers, one of us used Gravity’s Rainbow to cure some divergences appearing in Zero Point Energy 
(ZPE) calculations, at least to one loop[|]. The final ZPE result has been interpreted as an induced Cosmological 
Constant obtained as an eigenvalue of an appropriate Sturm-Liouville problem 1 . It is interesting to note that the 
same idea has been applied in a HL theory[l(J 2 in a FLRW background, the final result is that non trivial eigenvalues 
have been found depending on the parameters of the theory. Note that in GR, as we will show extensively in the 
next section, the cosmological constant cannot be considered as an eigenvalue of any Sturm-Liouville problem for the 
FLRW background in a mini-superspace approximation without matter fields. It appears therefore, that distortions 
of GR allow new results that otherwise should not be possible. It_remains to consider another distortion connected 
with the previous ones: a Varying Speed of Light (VSL) theory [l4| - |T9j . In this approach, one allows the speed of light 
to change in some specified way, in an attempt to solve the major cosmological issues of modern theoretical physics. 
It is well known, that one of the major features of Einstein’s theory of relativity is that the speed of light in a vacuum 
is always at constant rate. However, the cosmological problems that led to the theoretical introduction of dark matter 
and dark energy into modern cosmology have motivated some physicists to look for solutions in other directions, 
included the variation of the speed of light. In VSL, it is supposed that light travels faster in the early periods of 
the existence of the Universe and for this reason, it could solve problems related to the inflationary phase (flatness, 
horizon, homogeneity, etc....)fl7j }20l - l24j . Of course, this hypothesis breaks the Lorentz invariance. The VSL model 
has been embedded within the general framework of the time varying fine structure constant theory and reformulated 
as a dielectric vacuum theory [18j. Moreover, isotropy and homogeneity problems may find their appropriate solutions 
through this mechanism [25l-l27f. Recently quantum cosmological aspect of VSL models have been studied to see if 
the “Tunneling from Nothing” [28| and the “Hartle-Hawking No-boundary proposal” [29| can be better approached in 
this context©! [§l|. The purpose of this paper is to repeat the calculation of Ref. M in a VSL context to see if there 
are non trivial eigenvalues of an appropriate Sturm-Liouville problem, which will be interpreted as a Cosmological 
Constant induced by quantum fluctuations of the scale factor. The paper is organized as follows. In Sec. |H]we discuss 
the Wheeler-deWitt equations for Friedmann-Lemaitre-Robertson-Walker space time. While, in Sec. EU we show 
how it is possible to derive the Wheeler-deWitt equations for Friedmann-Lemaitre-Robertson-Walker space time in 
presence of varying speed of light. Conclusions are drawn in Sec lIVI 


1 See Ref. [gj for other applications in the context of Gravity’s Rainbow. 

2 See also Ref. EMI- See also Ref. fl3ll to see how Gravity’s Rainbos and HL theory can be connected. 
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A homogeneous, isotropic and closed universe is represented by the FLRW line element 

ds 2 = — N 2 dt 2 + a 2 (t) dfl 2 , (7) 

where 

dfl 2 = 7 ijdx l dxi ( 8 ) 

is the line element on the three-sphere, N is the lapse function and a(t) denotes the scale factor. Let us consider a 
very simple mini-superspace model described by the metric of Eq. ©. In this background, the Ricci curvature tensor 
and the scalar curvature read simply 


Rij ~ a 2 (t) 


7 ij 


and 


R = 


a 2 it) ’ 


respectively. The Einstein-Hilbert action in (3 + l)-dim is 

1 


S = 


167rG 


[ Cdt d 3 x = —[ Nyfg [K^Kij - K 2 + R - 2A] dt d 3 

JY.XI 167rG J SX 7 


(9) 


( 10 ) 


with A the cosmological constant, K t] the extrinsic curvature and K its trace. Using the line element, Eq. 0, the 
above written action, Eq. m, becomes 


S 



• 2 A s' 

a a — a H- a 


dt , 


where we have computed the volume associated to the three-sphere, namely V 3 = 2-7T 2 , and set N = 1. 
The canonical momentum reads 


( 11 ) 


<55 


37T 


^ = Td = -2G aa 


and the resulting Hamiltonian density is 


R = Tr a d — C 

G 2 37t 37t A 3 

= ~3^ 7r “ _ 4G a+ 4G 3° ' 

Following the canonical quantization prescription, we promote 7 r 0 to a momentum operator, setting 


da da 


( 12 ) 


(13) 


(14) 


where we have introduced a factor order ambiguity q. The generalization to k = 0, — 1 is straightforward. The WDW 
equation for such a metric is 


H'fy (a) 



d 2 _ qd_ 
da 2 a da 


Q ___ 

da 
97T 2 
4G 2 





T (a) = 0 . 


(15) 


It represents the quantum version of the invariance with respect to time reparametrization. If we define the following 
reference length a$ = -y/3/A, then Eq. (1151) assumes the familiar form of a one-dimensional Schrodinger equation for a 
particle moving in the potential 




(16) 
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with zero total energy. The potential U (a) resembles a potential well which is unbounded from below. When 
0 < a < ao, Eq. m implies U (a) > 0, which is the classically forbidden region, while for a > a 0 , one gets U (a) < 0, 
which is the classically allowed region. It is interesting to note that for for the special case of the operator ordering 
q = —1, one can determine exact solution]^. This can be easily verified by introducing the function 


1 9 ) — z o ( 1 9 


z(a) ={^G- 

where the solution can be written in terms of Airy functions, namely 

T (a) = aAi (z) + f3Bi (. z ). 

However, the wave function (ED. cannot be normalized in the following sense 

poo 

/ daa q ty* (a) $ (a). 


(17) 


(18) 


The same happens for the other special value q = 3. Even if the WDW equation (fT5l) has a zero energy eigenvalue, 
it also has a hidden structure. Indeed Ea. (fl5l) has the structure of a Sturm-Liouville eigenvalue problem with the 
cosmological constant as eigenvalue. We recall to the reader that a Sturm-Liouville differential equation is defined by 


d_ 

dx 

and the normalization is defined by 


p(x) 


dy (x) 
dx 


q (x) y ( x ) + A w ( x ) y (x) = 0 


/ dxw (x) y* (x) y (x). 
J a 

In the case of the FLRW model we have the following correspondence 

9<Z I 


(19) 


( 20 ) 


p{x) -> a q (t) , 


q (x) -> 


3n 

2G 

w (x) —> a q+i (■ t ) , 
y (x) ->• 'L (a) , 

A / 371- 

A W 2 G 


q -\-2 


(t), 


( 21 ) 


and the normalization becomes 


/ daa q+4x S>* (a) T (a). (22) 

Jo 

It is a standard procedure, to convert the Sturm-Liouville problem (|19[) into a variational problem of the form 3 

- fa dx v* i x ) [s (p (*) £) + <i (*)] y (*) 


F [y (x)] = 


f a dxw (x) y* (x) y (x) 


(24) 


3 Actually the standard variational procedure prefers the following form 

- \y* (x)p{x) -Ay (x)] a + fa dxp (x) (jj^y (x)) - q{x)y (x) 


F [y (x)] = 


fa dxw C 1 ) y* ( x ) y (*) 


(23) 


with appropriate boundary conditions. 
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with boundary condition to be specified. If y ( x ) is an eigenfunction of (IT9l) . then 

- la dx v* ( x ) l-t (p (*) -t) + 9 (*)] y ( x ) 


A = 


is the eigenvalue, otherwise 


Ai = min 

y{x) 


f a dxw (x) y* (x) y (x) 


~ I" dxy* (x) [£ (p jx) £) + g (a;)] y (x) 
J a dxw (x) y*(x)y (x) 


(25) 


(26) 


The minimum of the functional F [y (a;)] corresponds to a solution of the Sturm-Liouville problem (1191) with the 
eigenvalue A. In the mini-superspace approach with a FLRW background, one finds [l0| 


/ Vaa q 'k* (a) 


. Jl _ a d_ , 2 

9 ^ a da ^ 4 G? u 


* (°) 3Att 2 


f Vaa q + 4 \F* (a) ^ (a) 


4G 2 


(27) 


The best form of the trial wave function can be guessed by looking the asymptotic behavior of Eq. dTSl) . For a -» oo, 
we find 


d 2 q d 97 t 2 
da 2 a da 4G 2 


2 A z 
a - -a 


^ ( a ) — 


d 2 3Att 2 


da 2 4 G 2 


and when a —> 0, we find 

9 2 9tt 2 

Tfe 2 ~~ a 9a + 4G 2 


2 A , 
a - -a 


4-(a)~ 


d 2 _ qd_ 97^ 
Zk 2 “ + 4G 2 ' 


T (a) = 0 


T (a) = 0. 


(28) 


(29) 


When a — > 0, the previous equation can be exactly solved by a superposition of modified Bessel functions of the first 
(I v (a;)) and second kind (K u (x)). We findJH] 


to («) = 4 (gA) + C M <'-»>/ 2 A'„_ 1)/1 (ga 2 ) . 


(30) 


However, the solution (1301) is exact for a vanishing eigenvalue. Since our purpose is the evaluation of Eg. (1271) . we need 
a solution for a —> 0, which considers a generic not vanishing eigenvalue. This is described by 


T (a) = Ci exp - 


TG 2 exp - 


37ra^ 


8 G 


M 


q + 1 GE 2 q + 1 37T 


37ra" 

~8G 


U 


4 37t 2 4G 

<7 + 1 GE 2 g + l 37 t 
~4 3tT’ 2 ’4G ( 


(31) 


where M (a, 6, x) and U (a, b , x) are the Kurnmer functions. For practical purposes, it is useful to transform M (a, b, x) 
and U(a, b, x) in terms of I„ (x) and K„ (x). We find 


M(a + 1/2, 2 a + 1, 2a:) = T (1 + a) exp (x) I a (x) / {x/2)° 
U(a+ 1/2, 2 a + 1, 2a;) =exp (x) K a {x ) / (y/n (2a:) a ) 


(32) 


Since M(a,b,x ) is proportional to I a {x) which is divergent for large x, we will fix Gi = 0 to obtain normalizable 
solutions. Thus, we consider the following form 


’F (a) = exp ( — 


P a 2 


U 


<7 + 1 <7 + 1 


,P « 2 = 


(/3a 2 ) 


2 n(1-9)/4 




(Pa 2 

A (9—1)/4 ( — ]> 


(33) 


for the trial wave function and we plug (1531) into Eg. (1271) . After an integration over the scale factor a (t), one gets 


3A(/3)tt 2 _ f 0 °° dxxK 2 q _ 1)/4 (x) 


4 G 2 2/ 0 °° dxx 2 K 2 q _ 1)/4 (x) 


-B* + 


97T 2 

4G 2 ; 


(34) 
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where fi is a variational parameter and where we have rescaled the integrals with the help of the results of Appendix 
m By imposing that A (j3) be stationary against arbitrary variations of the parameter /3, we obtain 


This implies 


(B) = 3 fo° dxxK 2 (q _ 1)/4 (x) /2 G\ 2 d_ 
d/3 2 /“ dxx 2 I\ 2 q _ 1)/A (x) V 37r / d P 


Plugging (15S1) into Eq. dMl) . one finds 



A (/3±) 


4 


r (3+ 2) r(5_ £ ) 

r(^)r(V)' 


(35) 


(36) 


(37) 


It is easy to check that f}+ is a maximum and /3_ is a minimum. However /3_ is negative independently on <7 and this 
leads to a normalization (1221) in the range (— 00 , 0] which is non physical. A further exploration with a pure Gaussian 
choice, namely 


(a) = exp 



/3a 2 

— 


for q = 0 , leads to 


( n ) 

2 / Va-V* (a) 

_ Ji. + (3*A 2 a 2 

+ \2GJ U 

T(a) 

\2GJ 

f Vaa 4 'L* (a) T (a) 




(38) 


(39) 


The application of the variational procedure leads to imaginary solutions and therefore it will be discarded. It remains 
to test the following assumption 


W (a) = exp 



(40) 


suggested by the asymptotic behavior (1251) . In the next section we will discuss the choice CUl) as a particular case of 
a VSL theory. One could insist in this direction and try to explore other trial wave functions. However, choices (ED, 
EH]) and (|40l) have been chosen following the standard procedure for a variational approach. Therefore, the other 
choices can only be small variations of the proposed trial wave functions above mentioned. Therefore we are led to 
consider a distorted version of the gravitational field induced by a VSL theory. 


III. THE WHEELER-DEWITT EQUATION FOR THE 
FRIEDMANN-LEMAITRE-ROBERTSON-WALKER SPACE-TIME IN THE PRESENCE OF VARYING 

SPEED OF LIGHT 


A VSL cosmology model is described by the following line element 

ds 2 = — N 2 (t) c 2 (■ t ) dt 2 + a 2 (t) dfl 2 , (41) 

where dfl 2 is described by Eq.© and where c (t) is an arbitrary function of time with the dimensions of a [length/time]. 
The form of the background is such that the shift function N l vanishes. Thus, the extrinsic curvature reads 


SHj _ a jt) 

2 N(t)c(t)~ N (t) c (t) a (t) 9ij 


(42) 


where the dot denotes differentiation with respect to time t. The gravitational action fulfilling the Einstein’s Field 
equation with the speed of light explicitly written is 


S = T 


( 43 ) 
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where we have used the following relationship dx° = c(t) dt. It is easy to write the form of the reduced action of the 
mini-superspace. Indeed, reintroducing the speed of light into the action (flT)l) . one gets in (3 + 1) 

S = N(t) c 4 (t) ^ [K^Kij-K 2 + R- 2 A] dt d 3 x , (44) 


Using the line element, Eq. m , the above written action, Eq. dHD , becomes 


S = — 


37r 

4 G 



a 2 a — ac 2 ( t ) + ^a 3 c 2 (f) 
o 


dt 


(45) 


where we have computed the volume associated to the three-sphere, namely V 3 = 2i r 2 , and set N = 1. The canonical 
momentum reads 


- *JL 

5a 


37T 

2G 


aac 


it) 


(46) 


and the resulting Hamiltonian density is 


% = 7r a d — C 

G 


9 37T a 37T A o a . x 

W + iG 3“ c (i) 


According to the usual prescription where 7r a is promoted to an operator, we can write 


, , d 

7T a -IRC (t) — 


and introducing the factor ordering ambiguity 


tt; -> - (hc{t)Y a 

oa oa 


the WDW equation = 0 simply becomes 


Cj 9 37t a / \ 37r TV o /i / n t/\ /-x 

(t) + 4G3“ C (t) * (0) = a 


Following [IM3 , we assume that 


c ( t) = c 0 


a(t) 

ao 


(47) 


(48) 


(49) 


(50) 


(51) 


where ao is a reference length scale whose value will be fixed later. If the factor ordering is not distorted by the 
presence of a varying speed of light, one can further simplify the above equation to obtain 


d 2 _ qd_ 
da 2 a da 

where we have set N = 1 and the quantum potential is defined as 

A 


+ U c (a) T (a) = 0, 


U c (a) = 


3t r V 

2 Gh) 


a c b (t) 1 - —a = 


,_ 

\ 2 Ghal a ) 


= (_ 37rcg 


2+6a: 


A , 

1 - 


(52) 


(53) 












Note that the potential U c (a) vanishes in the same points where U (a) has its roots. Now, we are ready to discuss 
the analogue of Eg. (1271) in presence of a VSL distortion 4 . To this purpose Eg. (1571) can be cast into the form 


/ Vaa q T* (a) 


d 2 
da 2 


lA. 

a da 



\ 2 

( 3tt 

J a 2+6a 

\ 2lp a o a 


*(a) 


/ 'Daa q+4+6a 'i!* (a) d' (a) 


= 3A 


2lpd\ 


(56) 


where we have defined Ip = yjGH/ cq. Because of the VSL distortion, the asymptotic behavior of the trial wave 
function must be different compared to (127[) . Since 


K w (x) —ty/ it / (2x)exp (—x) x —^ oo 

Kq(x) -+ -ln(x) 

K v (x)->r(v)(x/2)~ v /2 


we find extremely useful the following assumption of the trial wave function 

, T . , \ -5+1 , a 3a / 

V (a) = a (pa) exp (-— 


(57) 


(58) 


which is a small variation of (|33|) . The exponential encodes the large a behavior, while (/3a) 3q encodes the small 
scale factor behavior and /3 is the variational parameter. Plugging (l58l) into Eo. (l56l) . after an integration over the 
scale factor a ( t ), we find 


where 


A,,„ (/?) = 3 


21 


2 3a 
P a 0 

37T 


3 ( 1 + 0 ) 


<?,(«)/3—+ 



(a) = i(g 2 - 24a - 2 q 7)T ■ 


We demand that 


dA 9jCl (/3) 
d/3 


3 /2^y 2 

2/3^ V 37r / 


^3 (1 + a) C q (a) /3 ^ 


37r 

21|ag“ 



= 0 , 


(59) 


(60) 


(61) 


where 


Pg (a) = - 


37T 


2ZpOQ a 


3(1 + a) Cg (ct) 


with the conditions 1 + a ^ 0 and 2 + 3a 7 ! 0. Plugging /3 into A giQ , (/3), one finds 

Ag,a (/3) = AF (a) ^ 2 1 + + 3 ; . 


(62) 


(63) 


Note that for the special case a = —2/3, one finds 


that it means 


#L-l±+( MM 2 fa- 2 - *(a) = 0, 

5a 2 a da ^2 Gha Q ~ J V 3 J J 


d 2 q d K 2 \ T , KK 2 T , , 

da 2 a da a 2 J 3 


(54) 


(55) 


where we have defined K = 37rag/ (2ip). This equation has exact solution in the form of a superposition of Bessel functions J„ ( x) and 
Y u ( x ). However to obtain eigenvalues one has to impose a large but finite boundary where the Bessel functions vanish. 
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The result is again dependent by the VSL parameter and on the reference scale a o. To this purpose we assume, 
without a loss of generality, that a 0 = kip. Then one gets 


-2 


Aq (&) — Ip 


37r 






2 T 3cy 
1 T Ot, 


,2 k 3a ) 3(1 + a)C q (a) 

To have one and only one solution, we find a stationary point for A q (a) and we impose that 

dAq (a) 


da 


= 0 


k = A(q, a) exp (B (q, a)) 


where 


and 


7T8 


= tt U 1 + a ) r “ 


3* 


1 T Sen. 


( -q 2 + 24a - 2q - 7) 


B (q, a) = 


(3 aq 2 - 72a 2 + 6 aq + 2 q 2 - 27a + 4 q+ 14) ^ + 48a + 32 


8 (q 2 - 24a + 2q + 7) 

In the table below, it is shown the result of the procedure (1551) for some specific choices of q 


q= i 

k 0 = 0.5779378002 

d = -2.007150679 

O 

II 

k 0 = 0.5843673484 

d = -1.988596177 

q = ~ 1 

k 0 = 0.6030705325 

d = -1.940190188 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


As depicted in Fig. (00), the couple (ko,a) does not represent the solution of the problem, because the point is stationary 



FIG. 1: Plot of Aq (a) as a function of a depicted for q = 1. The local minimum and the local maximum appear below the 
critical scale ko. For k = ko there is only a stationay point which disappears for k > ko. 

and not a local minimum. Rather we can interpret the couple (ko, ot) as a critical value below which a minimum and 
a maximum appear. In particular, as shown in Fig.. ©> for k < ko and a < a we have a minimum and for k < k 0 
and a > a, we have a maximum. In the spirit of the variational procedure only the minimum can be considered as 
the solution of the problem. Note that the lower the value of ko, the higher the value of A q (a). Note also that the 
value of ko is transplanckian. From the expression of A q (a), this is true when a > 0 or when a < —2/3, otherwise 
when —2/3 < a < 0, the behavior on k reverses. From the assumption (jSTjl . we see that for a > 0, we have c (t) Co 
when the scale factor a ( t ) ao and this is ruled out by observation. Therefore, the correct range of solutions is 
when a < —2/3. Furthermore to have also positive solutions we need a < — 1. The physical reason of why we obtain 
solutions in the negative range of a is that in the early universe one expects to measure strong quantum effects 
when the scale factor a ( t ) is really small. In this framework, this is realized with a speed of light which is really big 
compared to Cq- It is interesting to note that for k > ko, there is no solution at all and for k ko, A q (a) <C 1 in 
Planck’s units and a £ (d — e, a) with e > 0 arbitrarily small. 
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IV. CONCLUSIONS 


Motivated by the results obtained in different contexts of distorted gravity and in special way in the HL theorvjiol]. 
in this paper we have examined the possibility that the cosmological constant can be considered as an eigenvalue of 
an appropriate Sturm-Liouville problem of a VSL theory. This interpretation is not new and it has been explored 
in different contexts (8l-IToj. What is different in this paper is that the WDW equation on a FLRW background in a 
mini-superspace approach reveals a complete analogy with a Sturnr-Liouvillc problem and the cosmological constant 
has the natural interpretation of its related eigenvalue. The WDW equation has been examined taking account of 
the factor order ambiguity. We have probed ordinary gravity without matter fields with different families of trial 
wave functions and we have found no sign of a cosmological constant induced by quantum fluctuations. Of course we 
have not exhausted all the possible choices of the trial wave functions. However, the form we have adopted has been 
chosen using the standard criteria for a variational approach. Therefore, we can conjecture that for a mini-superspace 
approach without matter fields, a cosmological constant cannot be generated. The introduction of a VSL 

c<t>=c °(^r) (69) 


makes the situation different because of the power law on the scale factor. This modification is also supported by the 
following alternative definition of the speed of light 


c ( E/Epi) 


dE _ g 2 (E/Epi) 
dp C °gi (E/Epi)' 


(70) 


which can be easily extracted if one introduces Gravity’s Rainbow into the FLRW metric. In this formulation, the 
space-time geometry is described by the deformed metric 


ds 2 


N ' 2 (Jl dt 2 + a ' 2 W f/ o2 

gUE/E Pl ) dt + g^E/E Pl ) dU ^ 


(71) 


where g±(E/E P \) and g 2 (E / E P \) are functions of energy, which incorporate the deformation of the metric. Concerning 
the low-energy limit it is required to consider 


lim gi (E/Ep\) = 1 and lim g 2 ( E/E P \) = 1, (72) 

E / Epi~>0 E / Ep \—>0 


and thus to recover the usual background 0- Hence, E quantifies the energy scale at which quantum gravity effects 
become apparent. For instance, one of these effects would be that the graviton distorts the background metric as we 
approach the Planck scale. In a distorted FLRW metric the dispersion relation for a massless graviton is 


7^2 2 

E g i 


E 

~Ep 


2 2 

= P 92 


(£)• 


leading to m- Setting for example 


9i ( E/Ep \) — 1 
g 2 (E/E Pl ) = l+(^^j , 


(73) 


(74) 


one obtains a different, but equivalent form of the VSL. This formulation has the advantage to avoid technical 
complications as in Ref. 0- The choice in (l74l) appears to be connected also to the following potential 


a 4 ( t) 


6 


96nGb 34567 t 2 G 2 c 
a 4 (t) a 6 ( t ) 


(75) 


coming from a HL theory without detailed balanced condition. In this kind of potential, one discovers positive 
eigenvalues depending on the various coupling constants choices. However the potential (T75l) appears to be more 
flexible to produce positive eigenvalues. It is for this reason that the structure of the trial wave function we have used 
in this paper is more elaborated compared to a simple gaussian function which has bees used in a HL theory [To| . The 
procedure of finding a minimum for A (/?) of Eq. (1591) has produced a result depending on two other parameters, the 
power a and the reference scale k. A further minimization procedure allows to select one value compatible with the 
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procedure which however does not constitute the final answer, rather it has been interpreted as a critical value below 
which we have eigenvalues while above which we have none. Note that the appearance of eigenvalues compatible 
with the procedure is in the transplanckian regime and for negative values of a. Negative values of a have been 
found in Ref.jlij], even if the authors discuss the “Creation from Nothing” problem. Note that for Planckian and 
cisplanckian values of the scale ao, the eigenvalue does not appear and for larger scales, like the inflationary one, the 
whole expression in (l64l) becomes very small for every value of a < — 1. At this stage of calculation, we do not know 
if this behavior is simply a failure of the approach or further information can be extracted. 


Appendix A: Integrals for vH (a) = exp U( 




If the trial wave function assumes the form 


*(a)=exp l-^)U 


9 + 1 9 + 1 q _2 


4 ’ 2 


,Pa , 


then, for practical purposes, it can be transformed into 


2\ C 1 — ?)/ 4 


( P 0 ? 


, , (/3a 2 ) 

VSF K ^V\-T 

where we have used the identity (l32l) 

U(a + 1/2, 2 a + 1, 2x ) =exp (x) K a {x) / (tJtt (2a;)°) , 

with a = (q — l)/4 and x = /3a 2 / 2. Plugging the trial wave function (IA2I) into the kinetic term, one gets 

d 2 q d 


da 2 a da 


= —pa exp (-— U 


* («) 

/3a 5 


9+1 9+1 a 2 
—■ Sa 

If we multiply the expression (IA4I) on the left by a 9 '!'* (a) and we integrate over the scale factor, we find 

d 2 q d 


daa q 4-* (a) 


da 2 a da 


(a) = -P 2 J daa q+2 exp (-/3a 2 ) U 2 yy,/3a 2 ^ 


P~ 


daa 1)/41 2 


= --P 1 ^ [ dxxKf q _ 1)/4 (x) = -- -r 

\fpp2 a=s/x 71 JO n 


3 + q\ r f$-q 


where we have used the following relationship [37] 

/>00 

/ dxx a ~ 1 K^(cx)K jy (cx) 

Jo 


ia—3 


c a P (a) 


a + p, + is\ f a + p — v\ ( a + v — p\ ( a — p — v 


r 


r 


Rec > 0; Rea > |Rey| + |Re/r|. 


(Al) 

(A2) 

(A3) 


(A4) 


(A5) 


(A6) 


and where T (a;) is the gamma function. The contribution coming from the potential term without the VSL distortion 
is composed by 

J daa q 'i>* (a) a 2 ^ (a) = J daa q+2 exp (—/3a 2 ) \ Pa 2 ^ 


= 2 


P- 


yj (3/2a=y/x ^ Jo 


' , 3 ^2 ( 

daa K(g_x)/A ^ 


dxxKf q _ 1)/4 {x) = 

7rp 2 


3 + q\ ^ f5-q 


(A7) 
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and 


pOO pOO 

/ daa q ty* (a) a 4 4> (a) = / daa 9+4 exp (—/3a 2 ) £7 

Jo Jo 


2 , g +_1 £+1 0 2 

4 ’ 2 ,P 


4 


P~ 


^ r , s K 2 (p* 

~ Jo daaK (‘ 1 -^{ — 

dxx 2 K? q _ t)/i {x) = —LtT ( r f 7 ~ q 


y/p/la=v* n Jo w iy, ’~ ' 2/3 -t 1 V 4 / V 4 

The contribution coming from the potential term with the VSL distortion is composed by 

[ daa q <T (a) a 2+6a $ (o) = f daa q+2+6a exp (-/3a 2 ) C/ 2 (, /3a 2 
Jo Jo V 4 2 


and 


daa 3+6a Kl 


( Pa 2 


iq 1)74 V 2 ) y/pJla=V5n/3 2+3a+ 


ol+3 a r oo 

iET / dxx 1_ *" 3a K 2 q_ i) /4 (^) 

2 JO 


2 6 ° 


7t/3 2 5 +3«r (2 + 3a) 


r | 3 + 6a + r2 / 2 +3a A ^ (5 + 6a - q 


J daa q +* (a) a 4+6 “’I' (a) = J daa q+4+6ot exp (—/3a 2 ) U 2 ( g + , g + ,/3a 2 ) 




22+3( 


a /*oo 


y/j5/2a=y/x 7r(3' 


3+3a+- 


^ 2+3a ^(V 1)/4 ( a: ) 


22+6o 


TT^Sil+Sap (3 + 3a) 


p | 5 + 6a + q \ ^ ( 3 + 3a A T (7 + 6a-q 


(A8) 


(A9) 


(A10) 


Appendix B: Integrals for (a) = a 2 (/3a) 2 exp (— - 2 
If the trial wave function assumes the form 

T / \ _S±I , n / fiaJ 

T (a) = a 2 (/3 a ) 2 ex p f-- 

when we plug the trial wave function (IA2I) into the kinetic term, one gets 

A' 


a 2 _ 

9a 2 a 9a 
= 2/3a A (A 2 + p A — 2A) — a 2A /3 2 A 2 


* (a) = ^ exp (-^) ^a^ [( 9 2 - p 2 + 4p - 2q - 3) 


If we multiply the expression (IA4I) on the left by a 9 '?* (a) and we integrate over the scale factor, we find 


daa q +* (a) 


d 2 _ qd_ 
da 2 a da 


$( a ) = Ap + g2 ~ 2A ~ 2g + 1 r ( ~ 2 x +/ 7 


4A 


A 


(Bl) 


(B2) 


(B3) 


where T (a;) is the gamma function. The contribution coming from the potential term with the VSL distortion is 
composed by 


pOO pOO 

/ daa qx i>* (a) a 2+6a + (a) = / daa 1+6a+p P p exp (-/3a A ) 

Jo Jo 


1 \p— 6 cx—p —2 ( 2 + 6a + p 

=-+ > r 


(B4) 
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and 


n OO n OO 

/ daa q A>* (a) a 4+6a \l/ (a) = / daa 3+6a+p fi p exp (-/3a A ) 

Jo Jo 


A 


4 + 6 a + p 
A 


(B5) 
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